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Charge carrier dynamics in bulk M0S2 crystal studied by transient absorption 
microscopy 
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We report a transient absorption microscopy study of charge carrier dynamics in bulk M0S2 crystals at 
room temperature. Charge carriers are injected by interband absorption of a 555-nm pulse, and probed by 
measuring differential reflection of a time-delayed and spatially scanned 660-nm pulse. We find an intervalley 
transfer time of about 0.35 ps, an energy relaxation time of hot carriers on the order of 50 ps, and a carrier 
lifetime of 180 ± 20 ps. By monitoring the spatiotemporal dynamics of carriers, we obtained a diffusion 
coefficient of thermalized electrons of 4.2 ± 0.5 cm 2 /s, corresponding to a mobility of 170 ± 20 cm 2 /Vs. We 
also observed a time-varying diffusion coefficient of hot carriers. 
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I. INTRODUCTION 

Layered materials in which atomic sheets are stacked 
together by weak van der Waals forces can be used to 
fabricate two-dimensional systems. They represent a di- 
verse and rich, but largely unexplored, source of ma- 
terials. Graphene is the most well-known example*^— 
However, many other layered materials can also be ex- 
plored. Very recently, significant progress has been made 
on one type of layered transition metal dichalcogenide, 
M0S2. In 2010, two groups independently discovered an 
indirect-to-direct bandgap transition that occurs when 
varying the thickness from bulk to monolayer.^ In 2011, 
monolayer M0S2 transistors with a 10 s on/off ratio and 
a room-temperature mobility of more than 200 cm 2 /Vs 
have been demonstrated^ In 2012, several groups re- 
ported observation of valley selective optical excitation 
and luminescence in monolayer MoS2»^~— These studies 
paved ways to apply two dimensional M0S2 for photonics, 
electronics, and valley tronics. 

Studies on bulk M0S2 can provide complementary in- 
formation for understanding monolayers and their inter- 
action with environments and substrates, since the prop- 
erties of monolayer M0S2 are different from, but related 
to, the bulk. Here we report a transient absorption mi- 
croscopy study of charge carrier dynamics in bulk M0S2 
crystals. By time resolving the dynamics, we observed 
intervalley transfer, energy relaxation, and recombina- 
tion of carriers. Furthermore, spatial resolution of carrier 
dynamics allows us to observe diffusion of carriers, and 
measure their diffusion coefficients. These results provide 
fundamental information on charge carrier dynamics in 
bulk M0S2 crystals, and can be used to understand car- 
rier dynamics of M0S2 bulk and atomic layers. 
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II. EXPERIMENTAL SCHEME AND EXPECTED 
CARRIER DYNAMICS 



Figure Q] shows schematically the transient absorption 
microscopy setup. A passive mode-locked Tksapphirc 
laser is used to generate 100-fs pulses with a central 
wavelength of 780 nm and a repetition rate of 80 MHz. 
The output is used to pump an optical parametric os- 
cillator (OPO), which has a signal output of 1320 nm 
and an idler output of 1907 nm. We take a portion of 
the 780-nm Thsapphire output before entering the OPO, 
and combine it with the idler beam of OPO. By sending 
both beams to a beta barium borate (BBO) crystal, we 
generate their sum frequency with a central wavelength 
of 555-nm. This pulse is focused to the sample surface 
through a microscope objective lens with a spot size of 
1.6 fim (in full width at half maximum). It is used as 
the pump pulse. The probe pulse is obtained by sec- 
ond harmonic generation of the 1320 nm signal output of 
OPO in another BBO crystal, with a central wavelength 
of 660 nm. It is focused to the sample through the same 
objective lens, with a spot size of 1.3 /im. The reflected 
probe is collimated by the objective lens, and is sent to 
one photodiode of a balanced detector. A portion of the 
probe beam is taken before entering the sample, and is 
sent to the other photodiode of the balanced detector as 
the reference beam. 

The balanced detector outputs a voltage that is pro- 
portional to the difference between optical powers on the 
two photodiodes. In the measurements, we first block 
the pump pulse, and adjust the reference to match the 
reflected probe power, such that the balanced detector 
outputs a zero voltage. We then allow the pump pulse 
to reach the sample. It injects carriers, which change the 
reflection coefficient of the sample for the probe pulse. 
Hence, the balanced detector now outputs a voltage that 
is proportional to a differential reflection of the probe, 
Ai?/i?o- It is defined as the relative change of the probe 
reflection caused by the pump, (i? — Rq)/Rq, where R 
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FIG. 1. Schematics of the transient absorption microscopy 
system. 



and Rq are the reflection of the probe with the pump 
presence and without it, respectively. The advantage of 
using the balanced detector is to suppress the intensity 
noise of the probe pulse. Such a common mode noise is 
equally distributed on the two photodiodes, and hence is 
almost cancelled. To suppress other types of noise, we use 
a mechanical chopper to modulate the pump intensity at 
about 2 KHz, and use a lock-in amplifier to read the out- 
put of the balanced detector at that frequency. To mea- 
sure the differential reflection as a function of the probe 
delay, defined as the time delay of the probe pules with 
respect to the pump pulse, we change the length of the 
pump arm by moving the retroreflector in the pump arm. 
The time resolution in these measurements is about 250 
fs, since both pulses are about 180 fs at sample, mainly 
due to dispersion of the objective lens. To spatially re- 
solve the differential reflection signal, we scan the pump 
spot with respect to the probe spot at the sample sur- 
face by tilting the beamsplitter that reflects the pump 
into the objective lens. 

The M0S2 samples used for this study are natural 
occurring crystals from SPI Supplies. The background 
carrier density is below 10 15 /cm 3 . The samples were 
mounted on glass substrates. The surface of the samples 
is mechanically cleaved by using an adhesive tape in order 
to obtain a relative flat, clean, and fresh surface. All the 
measurements were performed under ambient conditions. 

The bandstructure of M0S2 is plotted schematically in 
Fig. [3J along with the expected carrier dynamics. The 
pump photon energy of 2.24 eV is large enough to ex- 
cite electrons near the direct bandgaps at K point, from 
both VI and V2 valanced bands to the CI conduction 
band. Once excited, the carriers can rapidly transfer to 
the lower valleys by intervalley scattering. These carriers 
then undergo an energy relaxation process, via emitting 
phonons, to reach thermal equilibrium with the lattice. 
Finally, the electrons and holes recombine. The dynam- 
ics is detected by using the 660-nm probe, with a pho- 
ton energy of 1.88 eV, within the A-exciton resonance^ 
Although the probe is resonant with the exciton transi- 
tion, we do not expect it detects real exciton population. 
Rather, the probe detects free carriers via carrier-induced 
nonlinearity at exciton transition^ In bulk M0S2, the 
exciton binding energy is about 25 meV, which is com- 
parable to thermal energy at room temperature.^ Hence, 
we do not expect a significant exciton population, and ig- 




FIG. 2. Bandstructure of M0S2 including the conduction (CI) 
and valence (VI and V2) bands. The vertical green and red 
arrows illustrate the pump and probe photon energies. The 
expected carrier dynamics is also shown, including interband 
excitation (1), intervalley transfer (2), energy relaxation (3), 
and interband recombination (4). 



nore their contributions. In fact, our experimental results 
also confirm this assumption (see Sec. IIVI) . 



III. DIFFERENTIAL REFLECTION TO PROBE 
CARRIER DYNAMICS 

In this section, we discuss in detail how to monitor the 
carrier dynamics with the probe pulse. For the normal 
incident probe pulse, the linear reflection coefficient of 
the unexcited sample can be obtained from Snell's law, 
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(n + l) 2 + , 



(1) 



where no and kq are the real and imaginary parts of the 
index of refraction. The latter, known as the extinction 
coefficient, is related to the absorption coefficient (ao) 
and light wavelength (A) by ko = (A/47r)ao- The charge 
carriers injected by the pump pulse can change the re- 
flection coefficient to 
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Since in our measurements the probe pulse is tuned to an 
excitonic resonance, the absorption is expected to play a 
dominant role. Hence, we assume real part is unchanged, 
n = no, and k = kq + Ak. Furthermore, if An only causes 
a small relative change of R, i.e. AR <C Ro, we have 
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In our measurements, AR/Rq is on the order of 10~ 3 . 
Hence, this condition is safely satisfied, and the differen- 
tial reflection 
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In our scheme, we probe the injected carriers by satu- 
ration of excitonic absorption induced by these carriers. 
In most cases, excitonic nonlinearities, such as absorp- 
tion saturation, are induced by real exciton populations 
via phase-space state filling. However, it has been estab- 
lished that excitonic transition strength can be changed 
by the phase-space state filling effect of free carriers, since 
the exciton wave function is composed of free-carrier 
statesJ^ Although this effects is generally weaker than 
those caused by exciton population, it is less selective 
to carrier energy distribution. Hence, it can effectively 
probe carriers with a high kinetic energy. Recently, we 
have shown that free carrier induced excitonic nonlin- 
earity can be used to probe ballistic transport of hot 
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and coherent plasma oscillations* in GaAs. 
In this effect, the change in extinction coefficient can be 
described by 
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where N is the carrier density, N s is the saturation den- 
sity and the cross section a is introduced to describe the 
effectiveness of carriers in saturating the probe transi- 
tions. Therefore, the differential reflection is related to 
the carrier density by 
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In the experiments, carriers are injected by interband 
absorption of the pump pulse. The on-axis peak energy 
fluence of the pump pulse, Fq, is related to its time av- 
erage power, P, by 
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where / is the repetition rate of the laser and w is the 
full width at half maximum of the focused pump spot. 
We assume that with each absorbed pump photon, one 
electron-hole pair is injected. The on-axis carrier density 
excited near the sample surface can then be written as 
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where ao and Tilo are the absorption coefficient and en- 
ergy of the pump photon, respectively. The lateral distri- 
bution of carriers is Gaussian, identical to the intensity 
profile of pump spot. From this procedure, we can con- 
vert the pump power to injected peak carrier density, by 
using the complex index of refraction of 5.4 + 3.1;/. 24 
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FIG. 3. (a) Differential reflection as a function of the probe 
delay (black curve and left axis). The inset shows the signal in 
early probe delays. The blue curve (right axis) is the carrier 
density deduced from the measured differential reflection, (b) 
Differential reflection as a function of injected carrier density 
for fixed probe delays of 0.2 (squares) and 2.0 ps (circles). 
The solid lines are calculated from Eq. [6] with the same value 
of N s = 3.7 x 10 19 /cm 3 . 



IV. RESULTS AND DISCUSSION 

The black curve in Fig. 0(a) (left axis) shows the 
measured differential reflection signal as a function of 
the probe delay, with a peak carrier density of 2.1 x 
10 19 /cm 3 . The inset shows the signal near zero probe 
delay. Clearly, the dynamics has a fast component of 
less than 1 ps and a slow one persists for several hun- 
dred ps. We attribute the slow process to recombination 
of carriers, i.e. carrier lifetime, and the fast process to 
the dynamical variation of a associated with evolution of 
carrier distribution in the conduction and valence bands. 
According to the bandstructure of M0S2, the pump pulse 
excites electrons from the valance bands to the K valley 
of the conduction band. Since the direct bandgap at K 
point is about 1.84 eV, electrons are excited with an ex- 
cess energy of about 100 - 200 meV. The pulse-width 
limited raising time of the differential reflection signal 
shows that these excited electrons change the excitonic 
transition strength immediately. This is an evidence that 
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exciton formation is negligible, since one would expect a 
relatively slow rise of signal associated with exciton for- 
mation. Since the Q valley in the conduction band is 
much lower than the K valley, these electrons transfer 
rapidly to the Q valley by emitting optical phonons. We 
attribute the initial decay of the signal in the first ps, 
with a time constant of about 0.35 ps, to this intervalley 
transfer process. This intervalley transfer time is compa- 
rable to electron phonon scattering time of other types 
of semiconductors (e.g. GaAs).— The Q-valley electrons 
still influences the exciton transition, as evidenced by the 
long-lived signal over 100 ps, but with a smaller cross 
section. Hence, we attribute the slow dynamics to re- 
combination of carriers in the Q valley. 

Since the carrier lifetime is on the order of 100 ps, we 
can assume that carrier density during the first a few ps 
is a constant, and equal to the injected density. We re- 
peat the measurement shown in the inset of Fig. |3]with 
various pump powers (and hence various injected peak 
carrier densities). The circles and squares in Fig. [3] show 
the differential reflection signal as a function of the in- 
jected carrier density measured at probe delays of 0.2 (the 
peak) and 2.0 ps, respectively. The dependence is non- 
linear, since the injected carrier density is rather high. 
Both data sets can be fit by using Eq. [6j with the same 
value of N s . By fitting many data sets with various probe 
delays, we deduce a value of N s = (3.7±0.1) x 10 19 /cm 3 , 
which is used to plot the two solid curves in Fig. [3jb) . 
With the known N s , we can deduce the carrier density 
as a function of probe delay from the measured differ- 
ential reflection. The result is plotted as the blue curve 
in Fig. [3Ka) (right axis). Hence, although the differential 
reflection is a nonlinear function of the carrier density, 
we can reliably deduce the latter from the former. 

Next, we spatially and temporally resolve the differ- 
ential reflection signal in order to study spatiotemporal 
dynamics of carriers. The carriers are injected by a fo- 
cused pump pulse with a Gaussian shape. Hence, the 
initial carrier distribution is 



N(r, 0) = iV exp[ 



(9) 



where <7o is the 1/e width of the pump spot. It is related 
to the full- width at half maxima, wq, by u>o = 2\/m2<7o- 
After injected, carriers diffusion in the sample plane 
and recombine. Although exciton formation is unlikely, 
due to the Coulomb attraction, the electron-hole pair 
still move as a unit in this ambipolar transport process. 
The dynamics is described by the ambipolar diffusion 
equation 
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dt t 

where D and r are ambipolar diffusion coefficient and 
carrier lifetime, respectively. Solving this equation, we 
have, 
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FIG. 4. (a) Carrier density as a function of probe delay and x 
over a long delay range of 500 ps. (b) A few examples of the 
spatial profiles. The probe delays are (from top to bottom) 
10, 30, 75, 100, and 200 ps, respectively, (c) Squared width of 
the spatial profile of the signal as a function of probe delay. 
The red line is a linear fit to the data points after 50 ps, 
which corresponds to a diffusion coefficient of 4.2 cm 2 /s. The 
blue line indicates the initial slope, which gives a diffusion 
coefficient of 18 cm 2 /s. (d) the carrier density measured at 
x = with the finite probe spot. The red line indicates a fit 
that gives a lifetime of 180 ps. 



where cr 2 {t) = a 2 +4Dt. Hence, we expect the profile re- 
mains Gaussian, with the squared width increase linearly 
with time. 

In order to quantitatively study the carrier dynamics, 
we measure the Gaussian spatial profiles along x axis for 
many probe delays, with a pump fluence of 10 /iJ/cm 2 . 
The results are plotted in Fig. Ufa). Clearly, the profile 
becomes lower and wider with increasing the probe delay. 
Figure Hfb) shows a few examples of the measured pro- 
files, with the corresponding Gaussian fits (solid lines). 
The squared widths deduced from the fits to all the mea- 
sured profiles are plotted in Fig. |U[c). After about 50 ps, 
the increase becomes linear, as expected from the am- 
bipolar diffusion model. From a linear fit (red line), we 
deduce a diffusion coefficient of 4.2 ± 0.4 cm 2 /s. We note 
that although the penetration depth is only on the order 
of 20 nm due to strong absorption, the diffusion of carri- 
ers along perpendicular direction is expected to be much 
slower due to the layered structure of M0S2. 

We are not aware of any previous studies on charge car- 
rier diffusion in bulk M0S2. However, from the measured 
ambipolar diffusion coefficient, we can estimate the car- 
rier mobility and compare with previous measurements 
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on mobility. Based on the fact that the effective masses 
of electrons and holes in M0S2 are similar^ we assuming 
that the mobilities of electrons and holes are compara- 
ble. Therefore, the ambipolar diffusion coefficient can be 
approximately treated as unipolar diffusion coefficients. 
Using the Einstein relation, D/ksT — fx/e, where ks, 
T, /i, and e are the Boltzmann constant, the tempera- 
ture, the mobility, and the absolute value of the electron 
charge, the measured D after 50 ps corresponds to a mo- 
bility of 170 ± 20 cm 2 /Vs. Early electrical measurements 
have deduced mobilities in the range of 100 - 260 cm 2 /Vs 
in nominally undoped bulk M0S2 crystals with carrier 
densities in the range of 10 15 to 10 16 /cm" Our result 
is reasonably consistent with these measurements. Fur- 
thermore, in the first 50 ps, the squared width increases 
sub-linearly, with a decreasing slope. The initial slope, 
indicated by the blue line in Fig. QJc), corresponds to a 
diffusion coefficient of 18 cm 2 /s. We attribute this to hot 
carrier effect. If the mobilities is not a strong function 
of temperature, hot carriers with larger thermal veloci- 
ties diffuse faster. Hence, the sub-linear increase of the 
squared width in the first 50 ps reflects the cooling of 
carriers, and provide a direct evidence of energy relax- 
ation of hot carriers. Such a hot carrier diffusion has 
been observed in other types of semiconductors, such as 

The spatiotemporal resolution of the carrier density 
also allow us to deduce the carrier lifetime. Figure |4](d) 
shows the measured carrier density as a function of probe 
delay with the probe spot overlapped with the pump spot 
(x = 0). The decrease of the carrier density is caused 
by both the carrier recombination and the transport of 
carriers out of the probing area. To deduce the carrier 
lifetime, we consider that the intensity of the Gaussian- 
shaped probe spot, with the center overlapped with the 
center of the carrier density profile, is 



J(r) = 7 exp[ -], 



(12) 



where Iq is on-axis peak intensity, and a p is the 1 /e width 
of the probe spot. By integrating the product of such a 
sensitivity function with the carrier density profile, we 
find that the carrier density measured by such a probe, 



N(x = 0) oc 
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Here, the first factor describes the effect of diffusion of 
carriers out of probing area, and the second factor ac- 
counts for recombination. We use Eq. [T3]to fit the data 
shown in Fig. QJd), with the known value of D. We ob- 
tain a satisfactory agreement with the data (red curve), 
and deduce a carrier lifetime of 180 ± 20 ps. We note 
that since M0S2 has an indirect bandgap, the radiative 
recombination of carriers is restricted by crystal momen- 
tum conservation. Accordingly, one would expect a long 
carrier lifetime due to radiative recombination. The short 
lifetime deduced here is probably limited by nonradiative 



recombination at defects and surface. It is also interest- 
ing to note that such a lifetime is consistent with that 
in atomically thin film M0S2 measured by time-resolved 
photoluminescence »2i 



V. SUMMARY 

We have used a transient absorption microscope to 
study charge carrier dynamics in bulk M0S2 crystals, and 
reveals several aspects of the charge carrier dynamics. 
We find that after the carriers are injected to the more 
energetic K valley, they transfer to the low energy Q val- 
ley with a time Constance of about 0.35 ps. After that, 
the carriers relax their energy in about 50 ps until they 
reach thermal equilibrium with the lattice. The rest of 
the carrier dynamics can be well described by ambipolar 
diffusion of thermalized carriers, with a diffusion coeffi- 
cient of 4.2 ± 0.4 cm 2 /s and a lifetime of 180 ± 20 ps. 
These results provide fundamental information on charge 
carrier dynamics in bulk M0S2 crystals, and can be used 
to understand carrier dynamics of M0S2 bulk and atomic 
layers. 
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